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Abstract. - We introduce stochastic driving in the Sznajd model of opinion spreading. 
This stochastic effect is meant to mimic a social temperature, so that agents can take random 
decisions with a varying probability. We show that a stochastic driving has a tremendous 
impact on the system dynamics as a whole by inducing an order- disorder nonequilibrium phase 
transition. Interestingly, under certain conditions, this stochastic dynamics can spontaneously 
lead to agents in the system who are analogous to Galam's contarians. 



Introduction. - The study of complex systems, in particular the application of statisti- 
cal physics methods to social phenomena, has recently attracted the attention of theoretical 
physicists [1-5]. From a statistical physics perspective, social systems are modeled as a col- 
lection of agents interacting through simple rules. In particular, the building (or the lack) of 
consensus in social systems has attracted much attention in recent years. A number of mod- 
els have been considered in order to mimic the dynamics of consensus in opinion formation, 
cultural dynamics, etc [6]. Among those models, Sznajd dynamics of opinion spreading has 
been subject of a great deal of work in recent years. Sznajd model is duly based on the trade 
union lema: "united we stand, divided we fall", and has been studied on different network 
topologies and for (slight) variations of the dynamics [6-11]. 

An important aspect of social and economic systems, recently discussed within opinion 
formation models, has been the presence of some agents called contrarians- namely, people 
who are in a "nonconformist opposition". That is, people who always adopt the opposite 
opinion to the majority [12-14]. In stock markets for instance, contrarians are those investors 
who buy shares of stock when most others are selling and sell when others are buying. The 
existence of a high proportion of contrarians in a society may play an important role in social 
dynamics (think for instance of referendums or stock market dynamics) [12]. 

In an attempt to include the contrarian effect in existing social models, a number of 
previous studies have considered contrarian agents as a initial condition, i.e. a given density 
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of contrarians is introduced in the model by hand for instance [12-14]. This is somewhat 
artificial and one would expect that simple models of opinion spreading should spontaneously 
lead to the existence of a fraction of contrarians among the population as some sort of emergent 
property. 

In this Letter we show that a contrarian effect can spontaneously emerge when stochastic 
driving is included in the model. As a typical model example we introduce stochastic dynamics 
on the Sznajd model. This randomness in the update of an agent opinion is meant to be 
a highly simplified description of the interplay between fashion/propaganda and a collective 
climate parameter, which is usually referred to as social temperature of the system [15-17]. We 
show here that social temperature in Snajd-type models leads to the spontaneous appearance 
of contrarian agents in the system. 

Mean-field approach. - Following the mean- field approach in [8] , we have considered the 
Sznajd model (the so-called "two over one" case in [8]), where two agents are chosen randomly 
and, if they are in consensus, then another randomly chosen agent is convinced by them. The 
Fokker-Planck equation (FPE) for the probability V(m, t) of having a "magnetization" to, at 
time t given a certain initial condition at time to < t is given by 



d_ 

dt 



V(m,t)= - ^(m(l-™>K*)) 



where N is the total number of agents. The magnetization density to = (N + — N~)/N 
measures the opinion state of the system and N + , N~ are the number of agents supporting 
the + or the — opinion, respectively (with N + + N~ = N). 

We now include a social temperature effect in the model by allowing the possibility that 
at every time step an agent follows the rules of the Sznajd model with probability p {jp < 1), 
while there is a probability 1 — p that those rules are not fulfilled (an agent adopts the opposite 
option than the one dictated by the rules). Then, for a given probability p, we arrive at a 
FPE for V p (m,t) that reads 

<9„ , , Id 



dt 



V p (m,t)= - — — ( y [(6p-5)m-(2p-l)m 3 )\v p (m,t)) 



+ G5v) ^([3-2 P -(2 P -l)TO^ p (TO,t)). (2) 
The stationary solution V^ at (m) results to be 

P?*( m ) « exp {^V jT ^ 2 ^g:jg du - ln[iV(3 - 2p - (2p - 1)to 2 )] } . (3) 

The analysis of this stationary solution for varying p shows that there is a threshold value, 
p = p c , such that for p > p c the system is bistable with a probability density Vf at {m) 
having two maxima at m± = ±v/(6p- 5)/(2p- 1). In this case the system gets ordered 
by spontaneously selecting one of the stable solutions m± . On the contrary, for p < p c the 
system becomes monostable and disordered with a magnetization density peaked at to = in 
which no dominant opinion survives. The threshold p c can be calculated in this mean-field 
approximation equating m + — to_ , at which all three extreme coalesce into a single minimum 
at to = 0, so that we find p c = 5/6. This behavior is shown in Fig. 1. 
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Fig. 1 - The effective mean-field potential obtained from Eq. For large values of p two minima 
m± appear. These move continuously toward each other as the threshold value p c = 5/6 is approached 
indicating a continuous phase transition. Bottom: Position of the potential minima vs p. When both 
extreme coalesce into m = 0, the transition point is reached. Top: Effective potential for (from left 
to right) p — 0.7, p = p c = 5/6, and p = 0.9. 



The picture emerging from the mean-field approach is clear. The effect of including thermal 
fluctuations in Snajd type models immediately leads to a contrarian-like effect. Some agents 
randomly take decisions that oppose the rules of the model, indicating some undecidedness 
in a fraction of the population. If such a fraction overcomes the critical threshold (1 — p c = 
1/6 ks 0.1667) the system will reach a stalemate situation, analogous to the contrarians effect 
discussed in [12-14]. 

Monte Carlo simulations. - In what follows we report on Monte Carlo simulations, in 
order to test the above discussed mean-field results. We have studied the model on regular 
lattices and small-world networks (which in the limit of high rewiring probability should 
reproduce the mean-field results). To make such an analysis, and in order to avoid the spurious 
antifcrromagnetic solution of the Sznajd original model, we have studied a convenient variation 
proposed in [18]: 

• rule 1': Chose an agent at random, say i, and if Sj x Sj +1 = 1, then Si-i and adopt 
the direction of the selected pair [i,i + 1], 

• rule 2': if Sj x Sj_|_i = — 1, then Sj adopts the direction of s,_x and Sj+x the direction of 

Si+2- 

Following [18], in case of disagreement of the pair (sj, Si+i), the rule 2' makes that the agent i 
"feels more comfortable" since it ends up with at least one neighbor having his same opinion. 
This variation of Sznajd model does not affect the basic behaviour and indeed has been shown 
to exhibit the same type of scaling features as the original model. 

We now introduce a stochastic mechanism in the dynamics as follows. At each Monte 
Carlo step we assume that, with a probability p, the rules are fully applied as indicated above, 
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Fig. 2 - Discontinuous phase transition for the one-dimensional lattice: the jump of the order param- 
eter m from ±1 to zero occurs abruptly. The discontinuity is more evident for larger lattices. Top: 
Stationary distribution of m for 8 = 0.44, 0.46, and 0.6 (from left to right), and system size N = 512. 
For small values of 8 a local maximum starts to develop at m = and becomes the global maximum 
for 9 > 8* . Data were averaged over 5 x 10 3 independent realizations. 



while the opposite option to the one dictated by the rules happens with a probability 1 — p. 
The probability p, in analogy with Weidlich [1,2] and Babinec [16], is defined according to 

p = Acxp ^ , (4) 

where a is some fixed parameter related to the strength of nearest neighbour interactions 
(we assume a > 0), and 9 is the collective climate parameter and plays the role of a (social) 
temperature. The normalization constant is A -1 = exp(a/0) + exp(— a/6). 
The asymptotic behaviour of p is: 

• if 9 — ► 0, we havep — > 1, indicating that without thermal fluctuations we recover Sanchez 
(and Sznajd) dynamics; 

• if 9 — > oo, we have p — ► 0.5, the probability of fulfilling the model rules or its opposite 
are the same. The model has a complete random behavior. 

Firstly we report on our results on the one-dimensional lattice. Each lattice site is occupied 
by one agent with opinion (spin) Sj G {+1,-1}. Starting from a random initial condition we let 
the system evolve towards its stationary state. For 9 — Q a consensus state arises (m± = ±1). 
However, as 9 is increased we observe transition towards the stalemate state as predicted by 
mean-field theory, nonetheless the transition is discontinuous (first order) as shown in Fig. 2. 

In order to compare with the mean-field results, we have studied our model in a fully 
connected network, which is expected to behave as a mean-field system. Indeed, in a fully 
connected network we observe that the transition between order and stalemate state becomes 
continuous. The qualitative agreement with mean-field results is apparent in Fig. 3. 
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Fig. 3 - Stationary distribution of m for a fully-connected lattice of N = 512 agents. Data were 
averaged over 100 independent runs. From left to right and top to bottom 9 = 0.30, 2.30, 2.60, and 
2.90. A second order transition towards a "stalemate" state is apparent: the most probably values of 
the order parameter m changes continuously from m = ±1 (bistable) to m = (mono stable) . 



We have also studied the intermediate cases between one-dimensional and mean-field (in- 
finite dimensional) limits by analyzing the model on small- world networks. Starting from a 
regular a one-dimensional lattice with periodic boundary conditions the links between neigh- 
bours are rewired with a certain probability r to a random site. Even for small values of r 
we observe a continuous phase transition from order to disorder as predicted by mean-field 
theory. We also observed that for increasing values of the rewiring probability r, the critical 
temperature, 8*, also increases, as shown in Fig. 4. 

The critical density of contrarians that it is required to reach the threshold (p c ~ 1 —p(0*)) 
is relatively large (the mean-field result p c ~ 0.33 in the case 2 against 2). However, such 
values are reasonable baring in mind that such a density of contrarians corresponds to a 
statistical average, dynamically generated by a large value of the social temperature (that 
is, each agent could sustain, convince or change its opinion dynamically). This mechanism 
is notably different from setting a fixed number or density p c of agents in a "nonconformist 
opposition" that will never follow the rules, as was done in previous studies [12-14]. In physical 
terms this difference is similar to the distinction between annealed and quenched disorder. 

Conclusions. - We have proposed a dynamical mechanism that leads to a contrarian-like 
effect analogous to the one described in [12]. However, in contrast to [12], we found that 
contrarians may spontaneously emerge from the dynamics when social temperature effects are 
taken into account. For low temperatures the system gets to a consensus where a majority 
opinion emerges just like in Sznajd type models. However, when temperature is above a critical 
threshold the density of contrarians is (on average) then high enough to make impossible for 
the system to reach a consensus and the opinion is equally divided between both options. 

Here, we have considered different forms and the most convenient prescriptions of the 
Sznajd model for the analytical and the numerical analysis. However, we have checked that 
the phenomenon is robust and does not depend on the particular form of the model. Moreover, 
since the Sznajd model (as well as many other two state opinion formation models) is similar 
to Ising type model [19] up to a certain extend, we can regard our results as a sophisticated 
manifestation in social systems of the ferromagnetic transition in spin systems. 
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Fig. 4 - The critical temperature 9* vs. the rewiring probability r. All simulations shown correspond 
to N — 512 agents and averages over 100 independent realizations. 

To conclude, the possibility of some external stochastic and/or deterministic influence on 
the agents of an opinion formation model, particularly regarding the possibility of some form 
of stochastic resonance [20], was recently analyzed by several authors [16,17]. It would be 
of great interest to study such an stochastic resonance effect (particularly its dependence on 
the size of the system [17]) in our model. This can actualy be done by including a fashion 
external field (for instance a periodic signal) combined with the noise effect coming from the 
social temperature. This is the subject of a forthcoming work. 

* * * 

We acknowledge financial support from Ministerio de Education y Ciencia (Spain) through 
Grant No. BFM2003-07749-C05-03 (Spain). MSL is supported by a FPU fellowship (Spain). 
HSW thanks the European Commission for the award of a Marie Curie Chair at Universidad 
de Cantabria (Spain). 

REFERENCES 

[1] Weidlich W., Sociodynamics-A systematic approach to mathematical modelling in social sciences 

(Taylor & Francis, London) 2002. 
[2] Weidlich W., Phys. Rep., 204 (1991) 1. 

[3] MOSS de Oliveira S., Oliveira P. CM. and Stauffer D., Evolution, money, war and com- 
puters (Teubner, Leipzig, Stuttgart) 1999. 

[4] Stauffer D., Introduction to Statistical Physics outside Physics - Physica A, 336 (2004) 1. 

[5] Galam S., Sociophysics: a personal testimony - Physica A, 336 (2004) 49. 

[6] Sec for instance Tessone C, Amengual P., Toral R., Wio H.S. and San Miguel M., Eur. 
Phys. J. B, 25 (2004) 403; Schneider J.J. and Hirtreiter Ch., Physica A, 353 (2005) 539, 
and references therein. 

[7] Sznajd-Weron K. and Sznajd J., Int. J. Mod. Phys. C, 11 (2000) 1157; Szna.id-Weron 
K. and Sznajd J., Int. J. Mod. Phys. C, 13 (2000) 115; Sznajd-Weron K., Phys. Rev. E, 66 
(2002) 046131. 

[8] Slanina F. and Lavicka H., Eur. Phys. J. B, 35 (2003) 279. 

[9] Stauffer D., Souza A.O. and Moss de Oliveira S., Int. J. Mod. Phys. C, 11 (2000) 1239; 
Ochrombel R., Int. J. Mod. Phys. C, 12 (2001) 1091; Stauffer D., Int. J. Mod. Phys. C, 13 
(2002) 315; Stauffer D. and Oliveira P. CM., Eur. Phys. J. B, 30 (2002) 587. 



Marta S. de La Lama, Juan M. Lopez and Horacio S. Wio: Spontaneous contrarian behaviour7 



[10] Stauffer D., Int. J. Mod. Phys. C, 13 (2002) 975. 

[11] Stauffer D., J. of Artificial Societies and Social Simulation, 5 (2001) 1; Stauffer D., AIP 
Conf. Proc, 690 (2003) 147; Stauffer D., Computing in Science and Engineering, 5 (2003) 
71. 

[12] Galam S., Physica A, 333 (2004) 453. 

[13] Stauffer D. and Martins J.S, Physica A, 334 (2004) 558. 
[14] Schneider J.J., Int. J. Mod. Phys. C, 15 (2004) 659. 

[15] Sanchez A.D., Lopez J.M., Rodriguez, M.A. Physical Review Letters, 88 (2002) 048701. 
[16] Babinec P., Phys. Lett. A, 225 (1997) 179; Kuperman M. and Zanette D., Eur. Phys. J. 
B, 26 (2002) 387. 

[17] Tessone C.J. and Toral R., Physica A, 351 (2005) 106. 

[18] Sanchez J.R., A Modified One- Dimensional Sznajd Model, cond-mat/0408518 (2004) . 
[19] Galam S., Europhys. Lett, 70 (2005) 705. 

[20] Gammaitoni L., Hanggi P., Jung P. and Marchesoni P., Rev. Mod. Phys., 70 (1998) 223. 



